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A New and Interesting Class of Limit Cycles 
in Recursive Digital Filters 

By V. B. LAWRENCE and K. V. MINA 
(Manuscript received February 1 , 1 978) 

Limit cycles oscillations often occur in recursive digital filters due 
to the quantization of products in the feedback section. A new and 
interesting class of limit cycles has been discovered and categorized 
for second-order sections that round either sign-magnitude or twos- 
complement products. These limit cycles are named rolling-pin limit 
cycles. They are completely defined by three integers and a simple 
construction rule and exist for B\ — B2 pairs lying within small 
rectangular regions in the B\ — B2 (coefficient) plane. Each set of 
integers completely defines the peak amplitude, the length, and the 
region of existence. The amplitude of these limit cycles can be made 
close to but does not exceed three times the Jackson peak estimate. 
Rolling-pin limit cycles often occur in filters with high Q poles located 
near dc or half the sampling frequency. When these large amplitude 
limit cycles occur, the idle channel performance of a filter is often 
unacceptable. Specialized techniques, requiring extra circuitry, can 
be used to suppress them. Alternatively, it may be less costly and 
more efficient to avoid the small rectangular regions within which 
the rolling-pin limit cycles exist in the B\ — B 2 coefficient plane. 

I. INTRODUCTION 

Oscillations often occur in recursive digital filters as a result of the 
nonlinear action of quantizing the products in the feedback sections. 
These oscillations occur in the least significant bits of the data and are 
called limit cycles.* These limit cycles influence the required internal 
data word length and hence the cost of the filter. 1 

This paper describes a new and important class of limit cycles that 
exist in second-order recursive digital filters. These limit cycles often 
occur in filters with high Q poles located near dc or half the sampling 



* These limit cycles are to be distinguished from the large limit cycles caused by 
overflow. 
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frequency. They are called rolling-pin limit cycles and derive this name 
from their characteristic shape when plotted in the successive value 
plane (Fig. 2a). The second-order section under consideration employs 
rounding of both feedback products in either sign-magnitude or twos- 
complement number format and is shown in Fig. 1. The rolling-pin 
limit cycles are defined by three integers, K, L, and M, and a simple 
construction rule. As seen in Fig. 2b, K is the constant step size in the 
handle of the rolling pin, L is the constant step size in the body of the 
rolling pin, and M is the number of steps of step size L. For each value 
of K, L, M, a unique set of limit cycles is completely defined. In this 
paper, emphasis is on the simplest case of K = 1. 

This class of rolling-pin limit cycles is important because of its 
unusually large amplitude. Unusually large amplitude limit cycles can 
lead to severe distressing tones in idle channel conditions. It will be 
shown that the peak amplitude approaches three times Jackson's peak 
estimate.* The concept of regions of the B x - B 2 plane within which 
complicated limit cycles exist is important. The existence of these 
isolated areas within which the various K, L, M limit cycles exist 
presents an interesting "patchwork quilt" look in the B\ — B 2 plane. A 
point of practical importance to note is that a small change in binary 
coefficient values (producing a pair of coefficients just outside the 
region of existence of rolling-pin limit cycles) can result in a 3:1 
reduction in ac limit cycle amplitude. By ac limit cycles, we mean limit 
cycles with period >2. 

Specialized techniques requiring extra circuitry 2 " 6 can be used to 
suppress rolling-pin limit cycles. These specialized techniques may 
increase the roundoff noise in the presence of a signal. Alternately, it 
may be less costly and more efficient to avoid rolling-pin limit cycles 
altogether, by avoiding the small rectangular regions within which 
they exist in the B\ — B 2 (coefficient) plane. 

In the sections that follow, we develop the explicit formulas for the 
peak amplitude, length, and regions of existence for the 1, L, M rolling- 
pin limit cycles. The various amplitude bounds and estimates that 
have been derived 7 " 12 in the past are examined relative to our exact 
peak values. A spectral analysis of rolling-pin limit cycles is also 
presented. A comparison of the roundoff noise and limit cycle power 
of a second-order section is made. 

II. PROPERTIES OF 1, L, M ROLLING-PIN LIMIT CYCLES 

It is important to get an overall view of the nature of this class of 
limit cycles before examining their detailed properties. For each value 
of L, M (the 1 will be omitted hereafter), a known number of rolling- 
pin limit cycles exist with identical bodies but handles that differ in 



* Jackson's estimate (Ref. 7) is the integer part of [0.5/(1 - B 2 )]. 
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Fig. 1— (a) Second order pole section, (b) Quantization characteristics. 



amplitude by one and overall length by four. For example, with L, M 
= 3, 4, there are three rolling-pin limit cycles with identical bodies but 
different peak amplitudes of 13, 14, and 15 and overall lengths of 34, 
38, and 42, respectively. The successive value plane plots (hereafter 
called the D\ — D 2 plots) for these three limit cycles are shown in 
Fig. 3. Each has its own cell of existence in the fii - B 2 plane. The 
three cells are horizontally contiguous as shown in Fig. 4. The bound- 
aries on B 2 are 



1 - (l/2)(l/9) > Bo > 1 - (1/2M1/8), 
and on B\ are 

2 - (3/2)(l/15) > \B X \ > 2 - (3/2X1/14) (amplitude 15) 



(1) 
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2 - (3/2) (1/14) > | fii| > 2 - (3/2) (1/13) (amplitude 14) 

2 - (3/2)(l/13) > \Bi\ > 2 - (3/2X1/12) (amplitude 13). (2) 

In Fig. 4, only negative values of B t (low-frequency poles) are shown; 
however, positive values of B\ also produce related 12 limit cycles. For 
convenience, we shall emphasize negative values of£i. The boundaries 
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defined in eqs. (1) and (2) are rational decimal (and rarely binary) 
numbers; however, in any implementation, the coefficients are binary. 
Accordingly, the equality signs in these constraints are usually not 
applicable. 
For other values of L, M, similar sets of nearly identical limit cycles 
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exist with either 3(L — 2) or 3(L — 3) contiguous cells of existence 
depending on whether the product L, M is even or odd. For each L 
there is a minimum value of M (but no maximum value) below which 
a type of degeneracy occurs in which the handle disappears into the 
body. Figure 5 shows the regions in the B\ — £2 plane for which each 
set exists for L = 3, 4, 5, 6 and four values of M for each L. In each 
case, the boundaries on B2 are of the form: 

1 - 1/[2(Q + 1)] > B 2 > 1 - 1/(2Q), (3) 

and the boundaries on B\ are 

2 - 3/(2P.) > |Bi| > 2 - 3/(2P,_,) (4) 

for i = 1, 2, • • • , (L — 2) or (L — 3). Both Q and P, are defined in terms 
of L and M in Section III. P, is the peak amplitude of the limit cycle. 
The allowed B\ and B2 values always correspond to complex conjugate 
poles. This can be proved using eqs. (3), (4), (6), and (12). 

Whenever either B2 or B\ is varied so that the outer boundary of a 
set of limit cycles is crossed, the rolling-pin limit cycle disappears and 
is replaced by a "normal" ac limit cycle with a smaller amplitude. The 
magnitude of the limit cycle is predicted by Jackson's estimate. 7 This 
effect should be borne in mind during the design of recursive filters, as 
an infinitestimal (binary) change in coefficient (and hence transfer 
function) can result in dramatic changes in the ac limit cycle properties. 
However, dc limit cycles larger than Jackson's estimate can still exist. 

This class of limit cycles has the unique feature that, given L, M, all 
states and properties are completely defined. In the following sections, 
the construction rules for rolling-pin limit cycles are given and formulas 
are derived for their peak amplitudes, length, and B\ — B2 boundaries 
in terms of L and M. 

III. CONSTRUCTION RULES FOR L, M ROLLING-PIN LIMIT CYCLES 

The construction rules for rolling-pin limit cycles are best under- 
stood by reference to Fig. 6, where a time sequence of an arbitrary 
rolling-pin limit cycle is shown. The 1, L, M rolling-pin limit cycles 
have amplitudes with half-period odd symmetry, i.e., 

Y\n + jj = -Y(n), (5) 

where N is the period (or length) of the limit cycle and is an even 
number. The rules differ slightly for even and odd values of the product 
L, M; thus they are discussed separately. 

3.1 Even L, M 

As evident in Fig. 6, a rolling-pin limit cycle consists of sections of 
constant slope (i.e., constant first difference) separated by smooth 
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Fig. 6 — Time sequence of 1, L, M limit cycle. 

transition regions (of constant second difference). Due to the odd 
symmetry, we need consider only half a cycle of the time sequence. 
Starting from the positive peak, we find that a zero crossing occurs 
during the decreasing section of constant step size L. This section has 
exactly M steps and is preceded and succeeded by smooth transitions 
to sections of constant unit step size. The transition section is of length 
L — 2. The handle and the body of the rolling pin are pictorially 
obvious in the D\ - D 2 plot in Fig. 2a; however, certain ambiguities 
arise in the discussion of the time sequence. First, every point in the 
D\ — Z>2 plot corresponds to a pair of values in the time sequence. 
Second, the transition region may be associated with either the handle 
or the body. The following convention will be used. The transition 
region is in the body of the rolling pin. The intersection of the handle 
and the transition region is also on the body. A value in the time 
sequence will be referred to as being on a certain part only when it is 
unambiguous. This would occur when both D\ — D2 points associated 
with the sample values were in the same part. 

For a given L, M, the bodies of all the constituent limit cycles are 
identical. It is the handle that varies in amplitude. Any single point on 
the body can be used to start the construction. We choose a value, Q, 
which is associated with the intersection of the handle and the body 
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(state A in Fig. 7). This value Q, labeled in the time sequence in Fig. 
6, is the next to last point on the unit step section immediately after 
the positive peak. For all even L, M, 



Q = LM/2 + L(L - l)/2 - 1. 



(6) 



The point Q defines the B2 bounds within which the limit cycle exists 
by 



1 >B 2 > 1- — . 

2(Q+1) 2Q- 



(7) 



Thus, for any even L, M we have defined the body of the rolling pin as 
well as the values of 5 2 for which it exists. Continuing on the decreasing 
portion of the limit cycle, the intersection of the body and the handle 
reoccurs at point B in Fig. 7. The A, D 2 values of B are —Q —2, — Q 
—1. To verify this, we start at A, with D\ = Q — 1 and counting to B 
using the construction rule just outlined, we have 



Df = Df - 2 - 3 (L-D-LM 



- (L - 1) - (L - 2) 2-1 



Df = Q - 1 ( -2 £ i ) + 1 - LM 




Fig. 7— Z>, - D 2 plot. 
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nR LM L(L-1) „ ___ r/T „ 
Df = -— +— — - — - - 2 - LM - L(L - 1) 
2 2 

Df--LM- L{L ~ 1) + l-2 

Z)f = -Q - 2. 

By symmetry, the states A' and 2J' also defined in Fig. 7 are — Q + 1, 
— Q and Q 4- 2, Q + 1 respectively. For the case of L, M = 3, 4 (the 
smallest set of rolling-pin limit cycles), Q = 8 and the bounds on B2 are 

1 - Vis > 5 2 > 1 - Vie. 

The common partial sequence in descending order is: 

Q 

12, 11, 10, 9, 8, 7, 5, 2, -1, -4, -7, -9, -10, -11, -12. 

state A state £ 

Having defined the body, we now define the allowed lengths of the 
handle (which fixes the peak amplitude of the limit cycle) and the 
bounds on B\. We do this by introducing the index J,. J, varies by unit 
steps from J max to Jmin inclusively, where 

Jnun = LM- L{L ~ 1] -3L + 5 (8) 

J max = LM- L(L " 1) -2. (9) 

The peak amplitude of the limit cycle is given by 

Pi t - Q + Ji. (10) 

The number of nearly identical rolling-pin limit cycles is: 

Jma* - Jmin + 1 = 3(L - 2). (11) 

Equation (11) has to be positive for a rolling-pin limit cycle to exist; 
hence, the minimum value of L is 3. The range on the peak amplitude 
Pi is between P™, and Pmin, where 

Pmax = Q + Jmax = 3 (^ - 1 J (12) 

•» min = *c ■ "min 

= 3 ^ - 3L + 4. (13) 

For L, M - 3, 4 
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5 < J, < 7 
13 < Pi < 15. 

The bounds on B\ for each of the L, M rolling-pin limit cycles are 
2--!">-£,>2--J- for i=l,2, ... ,3(L-2), (14) 

2.fi &ri-\ 

where 

Pi = P + i (15) 

and 

P = 3(LM/2 -L + l). (16) 

Crossing either the bounds on B 2 or the outer bounds on Bi results in 
elimination of the rolling-pin limit cycle. The amplitude of the resulting 
largest ac limit cycle will then be predicted by the Jackson estimate. 
The length of the limit cycle is N, where 

N = (AJi + 4L-6 + 2M). (17) 

Since N contains the term 4 J, and J, varies by 1 for each sucessive 
limit cycle in the set, N varies by 4. 

The example used throughout this section has been for L, M = 3, 4 
(the smallest of the set of the rolling-pin limit cycles). An example of 
a larger rolling-pin is L, M = 7, 8. This set has 3(L - 2) = 15 possible 
rolling-pin limit cycles. The peak amplitudes vary from 67 to 81. The 
values of Q and P Q are 48 and 66 respectively. The index J, varies from 
19 to 33. Figure 8 shows a magnified portion of the B\ — Bz plane 
where this set exits, and in Fig. 9 the time sequence of the largest and 
smallest constituent limit cycles are tabulated. 

3.2 Summary of construction rules and properties of the rolling-pin limit 
cycles 

3.2.1 Even L, M 

The results for even L, M rolling-pin limit cycles can be summarized 
as follows: 

(i) There are 3(L — 2) distinct limit cycles for any even L, M. Each 
has a time sequence with a half-period odd amplitude symmetry 

Y ( n+ lf) = _Y(n) - 

(ii) On the descending* half of the time sequence, there are 

(a) Ji + 1 steps of unit step size. 

(b) A smooth transition region (i.e., constant second differ- 
ence). 



* I.e., we start from the positive peak and scan the limit cycle. 
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Fig. 9 — Time sequences of L, M = 7, 8. 



(c) M steps of step size L. 

(d) A smooth transition region (also constant second differ- 
ence). 

(e) Ji — 1 steps of unit step size. 

(Hi) The peak amplitudes of the limit cycles are given by 

Pi = Ji + Q, 
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where 

Q = LM/2 + L(L - l)/2 - 1, 

J min = LM- L(L - l)/2 - 3(L - 2) - 1, 

J max = LM-L(L-l)/2-2. 

(iv) The region of existence of even L, M rolling-pin limit cycles in 
the B\ — Bi plane is defined by: 

2 - 75^- > -Sx > 2 - -J— for 1 = 1,2,... ,3(1,-2), 

Pi = Po + i, 
where 

Po = 3(LM/2 - (L - 2) - 1). 

3.2.2 Odd L, M 

Similar construction rules and bounds on B\ and B2 exist for odd 
values of L, M. These are simply summarized below. 

(i) There are 3(L - 3) distinct limit cycles for any odd L, M. Each 
has a half-amplitude period odd symmetry, 

Y(n + ^\ = -Y(n). 

(ii) On the descending half of the time sequence, there are 

(a) Ji + 1 steps of unit step size. 

(b) A smooth transition region (i.e., constant second differ- 
ence). 

(c) M steps of step size L. 

(d) A smooth transition region (i.e., constant second differ- 
ence). 

(e) Ji — 2 steps of unit step size. 

(Hi) The peak amplitude of the limit cycle is given by 



where 



Pi = Jt + Q, 



Q = (LM - l)/2 + L(L - l)/2 - 1 (18) 

Jnun = (LM - 1) - L(L ~ l) - 3(L - 3) - 1 (19) 



J™, = (LM-1)- L(L 2 l) -2. (20) 
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(iv) The region of existence of odd L, M rolling-pin limit cycles in 
the Bi — B2 plane is defined by 

1 "2(qTTT >B2 - 1 "^' 

• • , 3(L - 3), where 
Pi = Po + i, 

'(LM-1) 



for i = 1, 2, 



Po = 3 



- (L - 3) - 1 



(21) 



A point to note is that, for an odd L, M rolling-pin limit cycle to exist, 
we must have L > 5. 

The length of the limit cycle is N where 

N = (4J, + 4L + 2M-8). (22) 

Comparison of the equations for L, M odd with the equations for L, 
M even shows that the basic forms of the equations are identical if we 
interchange the following two quantities: 

LM LM-1 
2 2 

3(L - 2) <-► 3(L - 3). 

An example of odd L, M for L, M = 5, 7 is given. A magnified portion 
of the region of existence of this set in the Bi - B 2 plane is shown in 
Fig. 10. The time sequences of the largest and smallest members of 
this set are tabulated in Fig. 11. 



0.9815 



0.981 



0.9805 



0.980 

















»-!*» 


- 


48 


47 


46 


45 


44 


43 


-a a 4' 


-(2- 

1 


96' 


1 


1 


1 




-<2- 

1 



-1.968 



-1.964 



Fig. 10 — Region of existence of L, M = 5, 7. 
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IV. ADDITIONAL PROPERTIES OF 1,L,M LIMIT CYCLES 
4. 1 Peak amplitude comparisons 

In this section, the ratio, p, of the peak amplitude of a rolling-pin 
limit cycle and Jackson's peak amplitude estimate 7 are determined. In 
addition, the actual peak amplitudes of the L,M=3,4 and L, M = 3, 
12 rolling-pin limit cycles are compared to other calculated amplitude 
bounds or estimates. 7 " 11 

The peak amplitude of all rolling-pin limit cycles exceeds Jackson's 
estimate. It is now shown that the peak amplitude of a rolling-pin limit 
cycle approaches three times Jackson's estimate.* For convenience, 
we consider only even L, M. 

Jackson's estimate is 



2 (1-JB 2 ) 



(23) 



Since this depends only on B2, we can maximize p by minimizing JE. 
To do this, we pick the smallest value of B2 that will cause the rolling- 
pin limit cycle to occur. This value is given by 



so that 



B ° = 1 -2Q' 



2 2 



(24) 



(25) 



For a given L, M and a fixed B2 = 1 — 1/2Q, a set of rolling-pin limit 
cycles exist. To maximize p we pick the largest amplitude case [eq. 
(12)], 



* Parker and Hess conjectured in Ref. 8 that the limit cycle amplitude bound is three 
times the Jackson estimate. 
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\ z / 

Then 

LM-2 

For a fixed L, p may be made arbitrarily close to 3 by increasing M. 
For the two examples, L, M = 3, 4 and L,M=S, 12, p is 1.875 and 2.55 
respectively. 

Table I shows a comparison of the peak amplitudes for rolling-pin 
limit cycles of L, M = 3, 4 and L, M = 3, 12 with the bounds and 
estimates given by various authors. 7-11 As seen in this table, the actual 
peak amplitude is far less than the bounds predicted in Refs. 8 to 11. 
The Lyapunov bound 8 is the most pessimistic. The Sandberg-Kaiser 
peak amplitude estimate is V2/2 times their rms bound. 9 It is the best 
estimate for these examples. But even this bound is a factor of 6 higher 
for L, M = 3, 4 and a factor of 14 higher for L, M = 3, 12. 

4.2 Mean-square-value comparison 

Since all the states of the 1, L, M limit cycles are known, it is possible 
to develop an exact expression for the mean square value of the limit 
cycle sequence. The resulting expression is unfortunately complicated. 
The exact mean square value will however be compared to the Sand- 
berg-Kaiser mean square bound for a number of examples. We define 
the ratio y, 

y = Sandberg-Kaiser rms bound/actual rms value for various values 
of L and M. The Sandberg-Kaiser rms bound is 



SK = 



tf (i^k) >B " m 



(1 - B 2 ) Vl - B\/AB 2 

For Pi = Pmax, the above inequality on Bi and B 2 reduces to 

12Q > 4P + 3 or 2L(L - 1) > 1, (29) 

which always holds for rolling-pin limit cycles. Using the values of Bi, 
B 2 in the lower left-hand corner of the rectangle where the limit cycle 
exists, we have 

,„„ 2 _ 2(2P,) 2 (2Q) 2 (2Q-1) 

(&A) ~6(2P i )(2Q)-9Q-2(2P i r **" 

The ratio y was calculated for the largest limit cycle in the set 1, L, M 
for values of L = 4, 6, 8, 10, and 12 and values of M = Afmin to ( M mm 
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+ 10).* The results are plotted in Fig. 12. It can be seen that the SK 
rms bound is too high by a factor between 7 and 30. As either L or M 
is increased (the poles approach the unit circle), y increases approxi- 
mately linearly (pessimistically). 

4.3 Comparison of roundoff noise and limit cycle power 

In this section, the formula for the roundoff noise power from an 
isolated second-order section is derived in terms of P, and Q. A 
comparison between this noise power and the limit cycle power is 
given for L = 3, 4, 5, and 6 and four values of M for each L. It is 
concluded from the results of the comparison that the limit cycle 
power is likely to dominate. 

The average roundoff noise power at the output of the second-order 
filter section shown in Fig. 1 is given by 

<& = 2^(h(Z).W)y- < 31 > 



H(Z) is the transfer function of the second-order section. H(Z) is the 
complex conjugate of the transfer function. The quantization of prod- 
ucts in the presence of a signal is modeled by two white noise sources 
and is accounted for by the 2EI/12 term where E is the quantization 
step size of the filter. Evaluating the contour integral of eq. (31), we 
obtain 

* _ 9 El (1 + 1*2) 1 f32 . 

CTR - 2 l2(l-5 2 )(l + B 2 ) 2 -B?- (32) 

Expressing eq. (32) in terms of P, and Q, we have 

E 2 AP 2 O 2 

o% = 2^(4Q- 1) m _ 8Q) + 3Q 2 {8Pi _ 3 y (33) 

The limit cycle power is defined as 

N n? 

ol = E2 °^17' (34) 

where A are the values of the limit cycle. 
A useful comparison is to calculate 

/? = 10 log (ol/o%). (35) 

This gives an estimate of the relative importance of errors in the 
presence of a signal and errors in idle channel conditions. This is only 
an estimate, since the limit cycles have specific frequencies while the 
roundoff noise has a flat spectrum. In a cascade of sections, for 



* Afmin is the minimum value of M for which the rolling-pin limit cycles exist in their 

nondegenerate form. See Section V. 
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example, succeeding poles and zeros may completely eliminate or 
amplify one or more of the limit cycle frequencies. Spectral analysis of 
the rolling-pin limit cycles has been carried out. Figure 13 shows a 
typical result. These results indicate that several frequencies are 
present. Accordingly, the estimate /? is likely to be useful. In Table II, 
the value of /? in decibels is calculated for L = 3, 4, 5, and 6 and four 
values of M . In each case, the limit cycle noise power is larger. On the 
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Fig. 12 — Ratio of rms bounds. 
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Fig. 13 — Spectra of limit cycle. 
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Table II — Comparison of roundoff noise power and limit cycle power for 
various values of L, M 

Roundoff Noise limit Cycle fl-101og^dfl 

L, M Power o R Power ai ur 



3,4 


35.9 


108.6 


4.8 


3,6 


108 


264.8 


3.9 


3,8 


240.2 


488.6 


3.1 


3,10 


450.3 


780.1 


2.4 


4,4 


66.3 


228.24 


5.4 


4,5 


113.1 


365.5 


5.1 


4,6 


177.3 


533.5 


4.8 


4,7 


261.5 


732.4 


4.5 


5,5 


162.6 


582.4 


5.5 


5,6 


269.8 


915.8 


5.3 


5,7 


361.6 


1173.3 


5.1 


5,8 


533.9 


1624.9 


4.8 


6,5 


262.9 


993.4 


5.8 


6,6 


388.2 


1424.4 


5.6 


6,7 


545.9 


1924.8 


5.5 


6,8 


739.6 


2496.1 


5.3 



basis of these results, it is expected that limit cycle behavior will 
dominate in establishing the internal data word length of digital filters 
that generate rolling-pin limit cycles. 

4.4 DC and the small ac limit cycles 

A number of dc and small ac limit cycles are produced by the same 
second-order sections that produce the rolling-pin limit cycles. These 
small ac limit cycles are not accessible in the sense of Claasen et al. 12 
and their peak amplitudes do not exceed the Jackson estimate. The 
states and form of these small ac limit cycles can also be characterized 
by integers. Each of these small ac limit cycles has no handles. They 
only have a constant step-size body* with smooth transitional paths to 
their peaks. These smaller limit cycles can be constructed using the 
method outlined in Section V for rolling-pin limit cycles. 

The total number of these small ac limit cycles is ( Q — P D /3 — 1), 
where Q and P„ can be obtained from eqs. (6) and (16), respectively. 
The peak amplitude of the smallest ac limit cycle is ( P„/3 + 1). The 
largest amplitude of the small ac limit cycles is ( Q — 1). For L, M = 3, 
4 we have Q = 8 and P /3 = 4. Therefore, there are three small ac 
limit cycles. The smallest has a peak amplitude of 5 and the largest a 
peak amplitude of 7. 

The D\ — Z>2 plot of all the small ac and dc limit cycles for B\ = 
—1.875 and B 2 = 0.9375 (binary coefficients within the region for L, M 
= 3, 4) is shown in Fig. 14. The three large rolling-pin limit cycles are 



* For these various small limit cycles, the constant step size ranges from 1, 2- • • (L 
- 1), where L is the constant step size of the body of the associated rolling-pin limit 
cycle. 
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* 16 DC LIMIT CYCLES 
±1,±2,±3, ±4, ±9, 
±10, ±11, ±12, 
6 AC LIMIT CYCLES 

x x PEAK AMP 5 

c^— 3 PEAK AMP 6 
•■--• PEAK AMP 7 
a a PEAK AMP 13 




-- -15 
Fig. 14 — Successive value (D\ — Dz) plot small and large limit cycles. 



also shown in Fig. 14. The regions of existence of these small ac limit 
cycles together with the three rolling-pin limit cycles are shown in Fig. 
15. The rolling-pin limit cycles exist in a relatively smaller region. The 
small ac limit cycles exist in regions all the way up to the stability 
boundary B 2 = 1* We see from Fig. 15 that, as the rolling-pin 
boundaries are crossed, the larger rolling-pin limit cycles disappear. 
However, the small ac and dc limit cycles still remain. 

The number of small dc limit cycles can also be evaluated. This 
number is equal to 2P /3. Any state which satisfies the equation 



P P 



(36) 



is a dc state. These states represent the small dc limit cycles. The 
larger dc limit cycles can be evaluated from the following equation: 

Q + 1 < A = D 2 < Po. (37) 

For B, = -1.875 and B 2 - 0.9375, we have Q - 8 and P a = 12. 
Therefore, there are 8 small dc limit cycles, ±1, ±2, ±3, and ±4, all 



* Other limit cycles may also exist in the region outside the region of the rolling-pin 
limit cycles. No inference should be drawn on the existence of other limit cycles. 
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Fig. 15 — Region of existence of small and large limit cycles. 

satisfying eq. (36). In addition, there are 8 bigger dc limit cycles, ±9, 
±10, ±11, and ±12. These values can be obtained from eq. (37). These 
16 dc states are shown in Fig. 15. The state ±[Q, Q] = ±[8, 8], which 
should normally correspond to the Jackson estimate is not on any 
limit cycle. The next state after the state ±[8, 8] is on a bigger limit 
cycle with peak amplitude 13. This state ±[Q, Q] is unique because it 
is a state from which one can spiral out onto a bigger limit cycle. 

V. REGIONS OF EXISTENCE IN THE B, - B 2 PLANE 

Each member of a set of rolling-pin limit cycles exists in a rectangular 
region. The boundaries of the regions are rational numbers. As seen in 
Fig. 5, the members within a set exist in horizontally contiguous 
regions. In any real implementation, the coefficients are binary. 
Strictly, all statements refering to regions of existence should refer to 
the binary valued points within (and occasionally on the boundary of) 
the regions. For simplicity, we refer only to the continuous regions to 
demonstrate why the regions are rectangular; we develop the explicit 
boundaries for rolling-pin limit cycles. We now construct the values of 
Q, Ji, and P„ which were stated and used earlier. These variables, 
interestingly, depend on and follow from the assumed form of the body 
of the limit cycle which was constructed in Section III. The critical 
observation to make is that any state D h D 2 which generates the 
output Y [i.e., (Di,D 2 ) — > Y] defines a region in the B\ — B 2 plane. 
This region is an unbounded staircase of rectangles if D\D 2 ¥> 0, or is 
a semi-infinite slab if D\ or D 2 ¥* 0. The rectangles for the case of (3, 2) 
— » -1 or (-3, —2) —*■ 1 are shown in Fig. 16. The height and width of 
each rectangle is \1/D 2 \ and |1/A|, respectively. The center of each 
rectangle lies on the line defined by 
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Fig. 16 — Region of existence of a transition. Regions of existence include the solid, 
but not the dotted, boundary. 



B 2 D 2 +B1D1 Y, 



(38) 



with slope -D1/D2 and Si intercept —Y/D\. The B\ and B 2 intercepts 
of this line are both centers of rectangles. If D\ = 0, there is a horizontal 
slab of height 1 1/D 2 \ centered on B 2 = -Y/D 2 . For D 2 = 0, there is a 
vertical slab of width | 1/A| centered on Si = —Y/Di. This construc- 
tion was done for sign-magnitude products for the circuit shown in Fig. 
1. The rectangles include the solid boundaries but not the dotted 
boundaries since ±1/2 rounds to ±1. For twos-complement rounding, 
the same figure applies with minor changes in the dotted and solid 
lines since, in twos-complement, — 1/2 rounds to 0. Since any limit 
cycle is a sequence of transitions, the resulting region of existence is 
the region of the B\ — B 2 plane common to all the transitions. This 
region must be a rectangle. (When independent stability arguments 
are applied, the stability triangle is superimposed, which may further 
reduce the common region. This is the only condition that can cause 
other than a rectangular region.) 

Armed with this information, rolling-pin limit cycles were studied to 
determine the boundaries of their existence and to locate the critical 
transitions that set these boundaries. What was discovered was that 
the same transitions (relative to the overall shape of the successive 
value rolling-pin plot) always defined the B\ — B 2 regions. Small 
differences exist depending on the evenness and oddness of L and M. 
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For convenience, we consider M even. Since for sign-magnitude round- 
ing the states ( D h D 2 ) -* Y and (—A, -A) -» —Y result in the same 
B\ — B2 region, we need only consider the half of the limit cycle from 
the positive peak to the negative peak. For M even, the state —1 always 
occurs on the limit cycle. The body of the limit cycle is symmetrical 
with respect to this —1 point. The critical transitions that set the £2 
boundaries are shown in Fig. 17a. Two transitions ( T u and TJ) set the 
maximum value of £2 and one transition ( Tl) sets the minimum value 
of J3 2 . The D\ and D2 states in these transitions are defined in terms of 
Q = LM/2 + L(L - l)/2 - 1. The two transitions in the four 
consecutive states Q + 1, Q, Q — 1, Q — 3 completely determine the B 2 
boundaries. 

The boundaries on B\ are more complex since a horizontally contig- 
uous set of rectangles exist within each of which a separate rolling pin 
limit cycle occurs. The outer boundaries of this set of rectangles are 
determined by the two transitions T u and T L shown in Fig. 17b. The 
value of D\ in the transition that sets the upper limit on B\ is —(1 + 
ML/2 — L). The value of A in the transition that sets the lower limit 
on £1 is ML/2 — 1. The inner B\ boundaries that separate the nearly 
identical limit cycles in a set are fixed by the transitions at the peak of 
the limit cycle. Figure 17c illustrates this for a typical limit cycle in a 
set. As can be seen, the pair of consecutive transitions at either peak 
serves to set the inner boundaries on Bi. The peak transition of the 
largest limit cycle in the set sets the minimum bound on B\. This 



Q= LM/2+ L(L-1 )/2-1 





Q 3 Q-1 Q Q + 1 

---(LM/2-1) + (2L-3) 
(LM/2-1) + (L-1) 

(d) 

LM/2-1) 



-LM/2 A -LM/2 + I2L-1) 
-LM/2 + (L-1) 
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Fig. 17 — Critical transitions. 
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bound is identical to the ininimum bound on B\ determined by the 
body. In fact, this is how the peak limit cycle of the set is defined. The 
peak transition of the smallest limit cycle in the set sets the maximum 
bound on B\. This bound is identical to the maximum bound on B\ 
determined by the body. All peaks are possible between and including 
these two. The number of limit cycles in the set is determined in this 
way. 

When it is known which transitions set the boundaries and the 
values of the states in those transitions, the resulting boundaries are: 

1 > B 2 > 1 

2(Q + 1) 2Q 



2- 



2P, 



> - Bi > 2 - 



2P,--i 



for all i - 1, 2, 
and 



, 3(L — 2), where i indexes the consecutive states 



(LM 
Po = 3( — - L + 1 

Pi = P + i. 

Since Q and P, are integers, the boundaries are rational. The existence 
of the 2 in the denominator is expected since the boundaries must be 
associated with rounding products which are an integer ±1/2. The 
derivations of these boundaries is in principle possible by considering 
each transition in the body and locating a common region in the B\ 
— Bi plane. In practice, this is not accomplished; rather, the boundaries 
have been located using programming techniques. Insight into the 
process can be gained using the following approach. Assume for the 
critical Bi transitions in the body that B2 is effectively unity. (This 
must be rechecked later, but is in fact true.) The state equation for the 
circuit in Fig. 1 is then 



-D 2 + 



ft|2 -2X 



= Y. 



(39) 



The subscript R signifies sign-magnitude rounding. For the lower B x 
bound, the values of D u D 2 , Y are 

LM r n LM LM 

— + L-2. — -1, — (L + l). 



Using these values, eq. (39) becomes 



= LM - 3. 
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This nonlinear equation due to the rounding is valid for the interval 

3 1 (LM \ 1 

- > — ( 1 1 > - 

2 2X\ 2 / 2 

The lower Bi bound corresponds to the minimum algebraic value of 
[2 - (1/2X)] which occurs for the maximum value of X = (LM/2) — 1. 
This lower bound on B x is approachable but not attainable: 



B\ lower — — 2 + — " 



In a similar fashion, the values of A, D2, Y are 



2L - 1 - -=- I. I L - 1 - 



2 /' \~ 2 

This transition leads to the attainable bound: 



)• (-") 



Bi 



— 2 + - 



<¥-*«) 



Using a B\ between these limits and considering the two transitions in 
Fig. 17a, the bounds on B 2 are found to be: 

1 "2(QT1T >B2 - 1 "2Q- 

Checking back on the assumption that B 2 is effectively unity, for 

LM 
D 2 = _i___ 



1 \ / LM 



= -1- 



LM 



( LM\ 




(LM L(L - 1) 
2 ^ 2 + 2 


•■) 



ft' 



Since the niinimum L is 3, the R' operation yields a zero result and 
therefore B 2 is effectively 1. The R' operation is rounding except ±1/2 
rounds to 0. 

The quantity P, is the peak value of the ith limit cycle. The 
remaining unknowns are Jmin and e/ma X . They follow from Q and Pi. 

J min = Pi-Q = LM- U ~ -3 (L - 2) - 1 
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Jmax = P 3 (L - 2) ~ Q = LM - ^ ^ - 2. 

There are J max - Jmin + 1 = 3(L - 2) limit cycles in each set. 

This discussion can be extended to the cases "M-odd, L-even" and 
"M-odd, L-odd." The minor differences that occur arise from the 
following. For "M-odd, L-even," the value -L/2 - 1 is on the body of 
the limit cycle. This, however, leads to no changes in the equations 
just developed. For "M-odd, L-odd," the value -(L + l)/2 - 1 is on 
the limit cycle. This causes the minor variations in the equations 
summarized in Section 3.3. 

The final issue discussed in this section is a degeneracy that can 
occur for the smallest M for a given L. Essentially what happens is 
that the handle disappears into the body for some of the B\ rectangles 
nearest the origin. To determine if degeneracy occurs, the equation for 
e/ m in must be examined. For a given L, since there are «/, - 1 steps in 
the negative handle, if J min > 1, no degeneracy occurs. If Jmax > 1, but 
e/min is not, degeneracy occurs for that value of M. Figure 18 shows a 
plot of Jmin and Jmax as a function of M for L = 4 and L = 5. As can be 
seen, degeneracy occurs for L, M = 4, 3 and L,M=5, 4. 

VI. GENERALIZED K,L,M ROLLING-PIN LIMIT CYCLES 

Only the simplest class of rolling-pin limit cycles (1, L, M) has been 
analyzed in this paper. Other classes exist with values of K larger than 
1. The difference between the classes is K, which is the step size in the 
handle. The construction rules for the transition regions and the body 
are similar for all the classes. The other properties such as (i) regions 
of existence, (ii) peak amplitudes, (Hi) mean square value, (iv) number 




Fig. 18— Degeneracy criteria. 
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of cells, and (v) length can all be derived in terms of K, L, and M. The 
Di - D 2 plots, regions of existence, peak amplitudes, and number of 
cells for K, L, M =2, 6, 6 and K, L, M = 3, 8, 5 are shown in Figs. 19 
and 20, respectively. 

VII. CONCLUSIONS AND EXTENSIONS 

A unique set of unusually large limit cycles has been discovered and 
catalogued and are called rolling-pin limit cycles. The set exists for 
second-order feedback sections with two quantizers that round sign- 
magnitude or twos-complement products. The limit cycles are com- 
pletely defined by three integers K, L, and M, and a simple construction 
rule. The peak amplitude approaches three times Jackson's peak 
estimate. The limit cycles exist for Bi - B 2 pairs lying within rectan- 
gular regions in the Bi - B 2 plane. They occur often in filters with 
high Q poles near dc or half the sampling frequency. 

Specialized techniques requiring extra circuitry can be used to 
suppress rolling-pin limit cycles. These special techniques may also 
increase the roundoff noise in the presence of a signal. Alternately, 
these limit cycles may be avoided by making small changes in the 
binary coefficient values to produce a pair of coefficients just outside 
the region of existence of rolling-pin limit cycles and yet meet specified 
filter characteristics. 



60 



-K = 3. L = 8, M = 5 




60 



Fig. 19 — Successive value plot for generalized K, L, M limit cycles. 
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Fig. 20 — Region of existence of generalized K, L, M limit cycles. 

At present, these limit cycles are the largest (relative to Jackson's 
estimate) known in the region of the B\ — B 2 plane where |£i| > 1.875. 
For other regions, relatively large limit cycles have been found, but 
not systematized. It is expected that an approach similar to that 
presented will be useful. Other potentially useful extensions are for the 
cases of one quantizer and of truncation of products. 

At present, these results can be useful to digital filter designers, 
whenever poles with \B\ \ > 1.875 occur. Each design problem must be 
individually examined, however. A desirable goal is to incorporate 
these results into an automated design technique. 
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